Astronomy &: Astrophysics manuscript no. 0352 February 2, 2008 

(DOI: will be inserted by hand later) 



in 



o 



High-Resolution Near- Infrared Speckle Interferometry and 
Radiative Transfer Modeling of the OH/IR star OH 104.9+2.4* 

D.Riechers\ Y. Balega^, T. Driebe^ K.-H. Hofmann^ A. B. Men'shchikov^'^, G.Weigelt^ 

^ Mcix-Planck-Institut fiir Radioastronomie, Auf dem Hiigel 69, D-53121 Bonn, Germany 
^— N ' e-mail: riechers@mpifr-boim.mpg.de, driebe@mpifr-bonn.mpg.de, weigelt@mpifr-bonn.mpg.de 

("■^ ' ^ Special Astrophysical Observatory, Nizhnij Arkhyz, Zelenchuk district, Karachai-Chcrkessian Republic, Russia 

^\) , '' Institute for Computational Astrophysics, Saint Mary's University, Halifax, Canada 

^ , Received 27 February 2004 ; accepted 18 May 2004 

^: 

^ Abstract. We present near-infrared speckle interferometry of the OH/IR star OH 104.9-1-2.4 in the K band ob- 

^ . tained with the 6 m telescope of the Special Astrophysical Observatory (SAO). At a wavelength of A = 2.12 /im 

the diffraction-limited resolution of 74 mas was attained. The reconstructed visibility reveals a spherically sym- 

#^ , metric, circumstellar dust shell (CDS) surrounding the central star. The visibility function shows that the stellar 

^ ■ contribution to the total flux at A = 2.12 /im is less than ~ 50%, indicating a rather large optical depth of the 

fSJ ' CDS. The azimuthally averaged 1-dimensional Gaussian visibility fit yields a diameter of 47 ± 3 mas (FHWM), 

^\ , which corresponds to 112 ± 13 AU for an adopted distance of -D = 2.38 ± 0.24 kpc. To determine the structure 

^Z^ ' and the properties of the CDS of OH 104.9-1-2.4, radiative transfer calculations using the code DUSTY were per- 

>0 I formed to simultaneously model its visibility and the spectral energy distribution (SED) . We found that both the 

^— f^ ■ ISO spectrum and the visibility of OH 104.9-f 2.4 can be well reproduced by a radiative transfer model with an 

effective temperature Tcff — 2500 ± 500 K of the central source, a dust temperature Tin = 1000 ± 200 K at the inner 

shell boundary Rin — 9.1 i?* = 25.4 AU, an optical depth T2.2,jm ~ 6.5 ± 0.3, and dust grain radii ranging from 

r^ . fflmin = 0.005 ± 0.003 /im to Omax ~ 0.2 ± 0.02 /im with a power law n{a) oc a~^'^. It was found that even minor 

Oh' changes in Omax have a major impact on both the slope and the curvature of the visibility function, while the SED 

^ , shows only minor changes. Our detailed analysis demonstrates the potential of dust shell modeling constrained 

5_j ' by both the SED and visibilities. 

C^ . Key words, radiative transfer - techniques: image processing - stars: late-type - stars: AGB and post-AGB - 

IL" ' stars: OH/IR - stars: mass-loss - stars: circumstellar matter - infrared: stars - stars: individual: OH 104.94-2.4 

J^ 1. Introduction typically ~ Imag, a very small fraction varies irregularly 

■ - - ' with a low amplitude or does not show any detectable 

OH 104.9+2.4 (IRAS 22177+5936, AFGL 2885, NSV ^^^.^^ility. OH/IR stars are mostly low and intermediate 

25875) is an OH/IR type II-A class star. These objects ^^^^ (progenitor masses M < 9Mq), oxygen-rich single 

exhibit the maximum of their SED m the infrared (IR) ^^^^^ evolving along the upper part of the asymptotic 

around 6-10 ^m, while the 9.7 /im silicate feature is found ^-^^^ ^^^^^^ (^^B)^ qjj/jj^ ^^^^^ ^^^^^^ ^^^ sequence 

to be m absorption. In addition, they show strong OH ^^ ^^^-^^^ ^-^^^ variables towards longer periods, larger 

maser line emission at 1612 MHz/18 cm (Habing 1996) ^^^-^^^ ^^^^^^ ^^^ ^-^^^^^ mass-loss rates (Engels et 

and weaker emission in SiO (43.1GHz and 86.2 GHz), in ^^ ^ggg^ Repine et al. 1995). As a consequence of their 

H2O (22.2 GHz), and in the other hyperfine OH hues of ^.^^ ^^^ j^^^^ qjj/jj^ ^^^^^ ^^^ surrounded by massive, 

the same transition (1665/1667 MHz, 1720 MHz) (Herman ^^^.^^^^^ ^^^ geometrically thick circumstellar envelopes 

& Habing 1985). OH/IR stars are long-period variables composed of gas and dust. The strong maser emission 

(pulsation periods between 500 and 3000d) of variability ^^^^^^^ non-rotating CDS (Herman & Habing 1985), 

type Me, similar to long-period Mira stars. While the ^^-^^^ -^ ^^^^ ^^^^^^ ^^^^^^ obscures the underlying star. 

majority of OH/IR stars show a bolometric amplitude of ^, . ,, , -,11^, 

Observations as well as theoretical models of the struc- 

Send offprint requests to: D. Riechers ture, dynamics, and evolution of the atmosphere and CDS 

e-mail: riechers@mpifr-bonn.mpg.de of AGB stars have led to a detailed picture of these 

*Based on data collected at the 6 m BTA telescope of the objects. According to Habing (1996), OH/IR stars are 

Special Astrophysical Observatory in Russia characterized by the following parameters: a variation of 
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the bolometric flux Ajriboi > 0.7 mag, luminosities L > 
3 • IQ^ Lq, effective temperatures T^s < 3000 K, tempera- 
tures at the inner boundary of the CDS Tj^ — 1000 K, op- 
tical depths T9.7^in — 10^^ — 10^, the sihcate feature at this 
wavelength always in absorption, and high mass-loss rates 
(M ~ 10^^ — 10^** Mq/yt) at moderate outflow velocities 
{ve — lOkm/s). The AGB phase and therefore the phase 
of high mass loss typically lasts for some 10^ yr. Lorenz- 
Martins & de Araujo (1997) find T^s ~ 1800 - 2400 K, 
Tin — 650 — 1200 K, T9.7^m ~ 7 — 17, and the silicate fea- 
ture typically in absorption for many OH/IR stars. 

The mass loss is driven by both large-amplitude pul- 
sations and acceleration by radiation pressure. Therefore, 
correlations are found between the period and the terminal 
outflow velocity (e.g. Heske 1990). Detailed hydrodynam- 
ical models show that, due to the passage of shocks gen- 
erated by the stellar pulsation, the atmosphere is highly 
extended, thus enabling dust formation and the subse- 
quent acceleration of matter (see, e.g., Hofner et al. 2003, 
Schirrmacher et al. 2003, Jeong et al. 2004, and references 
therein). 

The assumption of a spherically symmetric dust shell is 
frequently motivated by the circularity of OH maser maps. 
High-spatial-resolution observations can give direct infor- 
mation about the extension and geometry of the CDS and 
therefore lead to more reliable constraints for the model- 
ing of the circumstellar environment, supplementing the 
information from the SED. Measurements of the near-IR 
(NIR) visibility may be used to determine the radius of 
the onset of dust formation as well as to constrain the 
grain sizes (Groenewegen 1997). 

OH 104.9-1-2.4 is a representative example of a type II- 
A OH/IR variable. From the variations of the 1612 MHz 
OH maser, a period of P = 1216 ± 54d was determined 
(Herman & Habing 1985), although values up to 1620 d 
have been found (Jones et al. 1990). The radial outflow 
velocity of v^ — 15km/s measured (te Lintel Hekkert 
et al. 1991) is in accordance with SiO and H2O maser 
observations (Engels et al. 1986). From maser phase lag 
measurements, an outer radius of the dust shell of i?s — 
(4.37 ± 0.42) • 10^^ cm was determined, while the absolute 
distance was found to be D = 2.38 ± 0.24 kpc (Herman & 
Habing 1985). A mass-loss rate of M = 5.58 • 10'^ Mg/yr 
was derived by Heske et al. (1990) . 

The paper is organized as follows: in Sect. 2, results of 
NIR speckle interferometric observations of OH 104.9+2.4 
are presented. Section 3 deals with the SED reconstruc- 
tion. The approach for the radiative transfer modeling 
(RTM), as well as the strategy of simultaneous mod- 
eling of SED and near-infrared visibility is outlined in 
Sect. 4.1. The method described was recently applied to 
other evolved stars, such as Red Rectangle (Men'shchikov 
et al. 1998), AFGL2290 (Ganger et al. 1999), IRC +10420 
(Blocker et al. 1999), NML Cyg (Blocker et al. 2001) and 
GIT 3 (Hofmann et al. 2001). Sections 4.2-4.4 give details 
of the modeling results, and Sect. 5 closes the discussion 
with a summary of the results, a conclusion and an out- 



look. Additional results will be presented in a subsequent 
publication (Riechers et al. 2005). 

2. Speckle interferometry observations and data 
reduction 

The if'-band speckle interferograms of OH 104.9+2.4 
were obtained with the Russian 6 m telescope of the 
Special Astrophysical Observatory (SAO) on September 
22, 2002. The data were recorded with our HAWAII 
speckle camera through interference filters with a cen- 
ter wavelength of 2.12 ^m and a bandwidth of 0.21 /xm 
(K' band). Additional speckle interferograms were taken 
for three unresolved reference stars (GJ 105.5, G 77-31, 
and HD 22686). With a pixel size of 27mas and a see- 
ing of 1.9 arcsec (1.5 • FWHM of centered long exposure), 
444+727object frames and 690+792 frames of the refer- 
ence stars were taken, each with an exposure time of 
268 ms. These interferograms were used to compensate for 
the speckle interferometry transfer function. The visibil- 
ity function of OH 104.9+2.4 was derived from the speckle 
interferograms using the speckle interferometry method 
(Labeyrie 1970). The reconstructed two-dimensional vis- 
ibility at 2.12 /im and the azimuthally averaged visibility 
profile is shown in Fig. 2. The azimuthally averaged vis- 
ibility decreases steadily to values below ~ 0.50 at the 
diffraction cut-off frequency (13.5cycles/arcsec). Thus, 
the CDS is almost fully resolved, and the contribution 
of the unresolved component to the monochromatic fiux 
at A = 2.12 /im is less than ^50%, indicating a rather 
large optical depth at this wavelength. In order to give 
a first, rough estimate for the apparent diameter of the 
dust shell, the azimuthally averaged visibility was fitted 
with a Gaussian and a uniform disk (UD) center-to-limb 
variation. We obtain a Gaussian fit FWHM diameter of 
47 ± 3 mas, which corresponds to a diameter of 112 ± 
13 AU for an adopted distance oi D ~ 2.38 ± 0.24 kpc 
(Herman & Habing 1985). The UD diameter is found to 
be 73 ± 5 mas, corresponding to 173 ± 22 AU. The quality 
of the visibility fits with both the Gaussian model and the 
UD are rather poor, as the shape of the observed visibility 
strongly deviates from that of a Gaussian or UD visibility. 
From our observations, no major deviation from spherical 
symmetry was detected, in accordance with the analysis 
of OH maser maps (Welty et al. 1987). Given the accu- 
racy and resolution of our visibility measurement, we find 
that the diameter ratio between major and minor axis is 
less than 8%. Therefore, we restricted the further discus- 
sion and the modeling itself to the azimuthally averaged 
1-dimcnsional visibility and the usage of a ID-radiative 
transfer code. 

3. Spectral energy distribution 

Photometric and spectroscopic data in the IR for 
OH 104.9+2.4 have been reported, for example, by 
Lebofsky et al. (1976), MerriU et al. (1976), Joyce et al. 
(1977), Ney et al. (1980), Price et al. (1983), Knacke 
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Fig. 1. Observations of 
OH 104.9+2.4 on September 
22, 2002, with the SAO 
6ni telescope. Left: 2- 
dimensional 2.12 /im 

visibility function of 
OH 104.9+2.4. Right: 

Azimuthally averaged vis- 
ibility of OH 104.9+2.4 in 
the K' band. 



et al. (1985), Herman ct al. (1986), Cobb et al. (1987), 
Jones et al. (1990), Le Squeren et al. (1992), Noguchi et 
al. (1993), Blommaert et al. (1993), Xiong et al. (1994), 
Jiang et al. (1997), Gonzalez- Alfonso et al. (1998), and 
Cutri et al. (2003). OH 104.9+2.4 was observed by IRAS 
(Infrared Astronomical Satellite) in 1983 and by ISO 
(Infrared Space Observatory) in 1996. In Fig. 3, the SED 
of OH 104.9+2.4 is shown. In addition to the cited data, 
recent data obtained by Yudin (2003) at the Crimean 
Astrophysical Observatory (CAO) on August 14, 2003, 
September 08, 2003, and December 10, 2003, in the J, 
H, K, L, and M bands (Tab. 1) are shown. The refer- 
ence star for these observations was BS8465 (J ~ 0.97™, 
H = 0.28"", K == 0.11", L = 0.00"", M = 0.20""). 
OH 104.9+2.4 shows strong long-term SED variations (see 
Fig. 3). The ISO data were most likely taken at mini- 
mum pulsation phase, as all other data points lie above 
the ISO spectrum. Therefore, we adopt $iso — 0.5. 
From the ISO data, a near-minimum bolometric flux of 
F^qP = 6.99 • lO^^^W/m^ can be derived. In order to 
find the bolometric flux at the phase of the SAO observa- 
tions, a bolometric flux at maximum phase has to be de- 
termined first. To obtain this bolometric flux, we defined a 
maximum envelope representing the maximum-luminosity 
SED as depicted by the black line in Fig. 3, which takes 
into account the wavelength dependence of the pulsation 
amplitude, and obtain i^^f = 2.30 • 10"^" W/m^. With 
K-h&nd data measured by Jones et al. (1990), Noguchi 
et al. (1993), Cutri et al. (2003), and Yudin (2003) and a 
simple cosine-like pulsation model, we derived a pulsation 
period of P = 1500 ± 11 d, applying a 4-parameter least- 
squares fit (/(x) = a • cos{2ttx / P — b) + c), as shown in 
Fig. 3. We would like to point out, that this model is not 
physically accurate, but a suitable choice regarding the 
low number of data points. With this model fit, we find 
$SAO - 0. Therefore, F^^^ ^ 2.30 • 10^^° W/m^ can be 
assumed, corresponding to a factor of ^ 3.3 relative to the 
flux at $iso- This difference in Fboi has to be taken into 
account in the modeling process. 



OH 104.9+2.4 is highly reddened by the interstellar 
medium (ISM) and its CDS. The interstellar reddening 
was taken into account by adopting the method of Savage 
& Mathis (1979) with Ay =■ 3.1 • E{B - V). For a dis- 
tance of D = 2.38 ± 0.24 kpc, we found Ay = 4.09. 
This correction was applied to all models as described by 
Men'shchikov et al. (2002a). 

4. Dust shell models 

4.1. The radiative transfer code 

The radiative transfer calculations were performed using 
the pubhcly available code DUSTY (Ivezic et al. 1999). 
The main assumptions involved are (for single-shell mod- 
els): (i) constant spherical symmetry with p{r) oc r~", (ii) 
a constant dust-to-gas ratio, (in) a constant dust opacity 
and size distribution, (iv) radiative equilibrium, (v) the 
same temperature rd(a,r) = Td(r) for all dust grains, 
and (vi) isotropic scattering. With these assumptions, the 
only relevant property of the input radiation is its spectral 
shape f\ = F\/F. The overall optical depth of the dust 
envelope at some reference wavelength t\ fixes the solu- 
tion, allowing only for normalized density distributions pd 
describing the spatial variation of the dust, and the wave- 
length dependence of the optical properties of the dust 
grains. 

The DUSTY formulation of the radiative transfer 
problem for a dusty envelope is very suitable for mod- 
cling of IR observations, since it minimizes the number of 
independent model input parameters. The only parame- 
ters that have to be specified are (i) the spectral shape 
of the central source of radiation f\ (i.e. the variation 
of the monochromatic fiux with wavelength), {ii) the ab- 
sorption/emission and scattering efficiencies of the grains 
(i.e. the chemical composition and the grain-size distribu- 
tion), (Hi) the normalized density distribution of the dust 
p{r)/pin, (iv) the relative thickness of the dust shell (i.e. 
the ratio of outer to inner dust shell radius, Tout/nn), (v) 
the dust temperature Tin at the inner boundary rjn, and 
(vi) the overall optical depth at a reference wavelength t\. 
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Fig. 2. SED of OH 104.9+2.4. The ISO 
SED appears to be taken at near-minimum 
pulsation phase ("I>iso = 0.5). The black 
line depicts a maximum envelope for the 
SED, which was used to derive the bolo- 
metric flux at maximum phase. 
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Fig. 3. Recent observations of OH 104.9+2.4 at A = 
2.2 /im. Measurements by Jones et al. (1990), Noguchi 
et al. (1993), Cutri et al. (2003) and Yudin (2003) are 
shown. A 4-parameter cosine function is fitted (least- 
squares method) to the data and provides a pulsation 
period of P = 1500 ± lid. The left vertical line in- 
dicates the date of the ISO measurements and sub- 
stantiates the assumption of $iso ~ 0.5. The right 
vertical line indicates the date of the SAO measure- 
nents and allows us to derive $sao — 0. 



Table 1. Infrared photometry of OH 104.9+2.4. The observations were carried out by Yudin (2003). All observations are given 
in magnitudes. 



Date 


Julian Date 


J 


H 


K 


L 


M 


Aug 14, 2003 


2452866.50 




9.04 


5.20 


1.16 


-0.06 


Sep 08, 2003 


2452891.42 


11.53 


9.19 


5.32 


1.24 


-0.06 


Dec 10, 2003 


2452984.20 






5.63 


1.54 


0.18 



For a given set of input parameters, DUSTY iteratively 
determines the radiation field and the dust temperature 
distribution (generalized output quantities) by solving the 
equation of radiative transfer. 



4.2. Previous modeling of OH 104.9+2.4 

There are two previous studies which deal with dust shell 
modehng of OH 104.9+2.4. The first dust sheU models 
were reported by Ivezic et al. (1995) using a previous 
version of the DUSTY code. It was mainly based on 
IRAS data. The authors found Fboi = 7.6 • 10"" W/m^, 



L = 

T'2.2 fim 



1.19 • 10"* Lq for D = 2.30 kpc, v^ = 15.1km/s, 
= 1.12, and M = 3.87 • 10"^ Mq/jt. 



A recent two-component dust shell model was reported 
by Suh (2002), using only the ISO SWS/LWS spectra as 
observational constraints. The author introduced a com- 
bined envelope/disk model to improve his fit to the ob- 
servational data. The reference optical depth for his best- 



fitting model was tiq 



/ini 



10.0. To date, there has not 



been any modeling of OH 104.9+2.4 that uses observa- 
tional constraints from both the SED and visibility mea- 
surements. 
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For the analysis of OH 104.9+2.4 based on the radiative 
transfer modeling, ^ 10^ models (discussed in Sect. 4.4) 
were calculated. To simultaneously fit the SED and vis- 
ibility data (i.e. with identical model parameters), -Fboi 
had to be adjusted for the epoch of the visibility mea- 
surements, since the ISO spectrum and our visibility data 
were taken at different pulsation phases of OH 104.9+2.4 
(see Fig. 3). Therefore, the best- fitting model (Figs. 4 and 
5) consists of an SED model assuming the 1996 ISO flux 
and a visibility model assuming the 2002 SAO flux. The 
results are all summarized in Tab. 2. 



4.3.1. SED fit 

The final model that fits the ISO SED data of 
OH 104.9+2.4 best is presented in Fig. 4 (left panel). 
The black-body (BB) effective temperature of the cen- 
tral source of radiation is Toff = 2500 K. In addition. 
Tin - 1000 K, pir)/pin ex r'^-" and rout/nn - 10^. 
Assuming the overall influence of the dust formation zone 
on the SED to be rather small, a standard n{a) oc a~^'^ 
grain-size distribution (GSD) function (Mathis, Rumpl, 
& Nordsieck 1977, hereafter MRN) was chosen, with a 
ranging between amin = 0.005 /nm and amax = 0.2 /im. 
A dust composition with 95 % of the optical constants 
for warm silicates from Ossenkopf, Henning, & Mathis 
(1992, hereafter OHM) and 5 % astronomical silicates from 
Draine & Lee (1984, hereafter DL) was found to give the 
best results. The OHM warm, oxygen-deficient silicates 
are based on observational determinations of the opacities 
of circumstellar silicates as well as laboratory data and 
also take into account the effects of inclusions (amorphous 
and crystalline substructures as well as other materials) 
on the complex refractive index, especially at A — 8 /im. 
The reference optical depth in our best-fitting model is 
T2.2fim = 6.5, corresponding to rg.r^^m — 14.0 for this 
model at the prominent silicate absorption feature. The 
observational value of i^boi = 6.99 • 10~^^ W/m^ from the 
ISO data can be reproduced well by our model. 

From the shortest wavelength at A = 2.38 /im up to 
A = 196 /im, the model provides a good fit to the obser- 
vations. The location, shape, and strength of the silicate 
feature at A = 9.7 /im caused by SiO stretching vibrations 
is fairly well reproduced by the given composition. In the 
A ~ 18 /im region, where Si02 bending vibrations dom- 
inate the spectrum, there is a noticeable deviation from 
the observational data. Above A ~ 100 /im, the model falls 
off too steeply to fit the data. On the one hand, this may 
be caused by deteriorating data quality towards the long 
wavelengths. On the other hand, the discrepancy at longer 
wavelengths may be due to the fact that DUSTY uses a 
GSD as input but internally averages dust properties over 
the entire size distribution. 



From Fig. 3 it is obvious that the ISO SED and the SAO 
visibility data were taken at different pulsation phases. 
While the ISO observations were made at near-minimum 
phase, the SAO data were taken at near-maximum phase. 
Therefore, it is quite natural that the two constraints can- 
not be fulfilled by one model alone with the same Fboi- A 
factor of ^ 3.3 in terms of the flux at ISO phase has to be 
introduced to the visibility model if it is to be comparable 
to the SED model, as outlined in Sect. 3. If we combine 
the best-fitting SED model with i^boi = 6.99 • 10"" W/m^ 
from the previous subsection with a corresponding visibil- 
ity model with Fboi = 2.30- 10"^" W/m^, we obtain a very 
good simultaneous model for the SED and the visibility 
of OH 104.9+2.4, as shown in Fig. 4. However, it should 
be stated that the physical conclusions that can be made 
based on this simultaneity itself are limited because of the 
simplifications involved in the phase alignment and the 
simplicity of the pulsation model itself (only flux scaling, 
assumption of time-independent Toff and i?*). Therefore, 
all futher derived quantities are given for the phase of the 
ISO measurements unless explicitly stated otherwise. 

4.3.3. Further model results 

The top left panel of Fig. 5 shows the fractional flux con- 
tributions of the emerging stellar radiation, the scattered 
radiation, and thermal dust emission as a function of wave- 
length, for our best-fitting model. At 2.2 /im, the fiux is 
dominated by thermal dust emission (45 %) and scattering 
(35%). Direct stellar light (attenuated flux) contributes 
only 20% at this wavelength. Thus, scattering contributes 
nearly to the same extent to the radiation of the CDS in 
the K band as thermal dust emission. At shorter wave- 
lengths scattering dominates, while dust emission is the 
major flux contribution at wavelengths longer than 3 /im. 
The attenuated fiux from the central source plays a mi- 
nor role at all wavelengths under consideration (0.01 /im - 
0.03 m) due to nearly total obscuration. 

The top right panel of Fig. 5 shows the normalized 
intensity distribution at 2.13 /im as a function of an- 
gular distance. The barely resolved central peak corre- 
sponds to the central star. At the given optical depth of 
T2.2 /^m — 6.5, there is no evidence for limb-brightening 
effects at the inner boundary, due to the extreme opti- 
cal thickness of the CDS. If we define an outer radius of 
the dust shell i?out° ^^ the point where the intensity at 
2.13 /im is /out~^" — 10"^° • /max, this point is found at 
approximately i?out^ — 5.8arcsec. 

The bottom left panel of Fig. 5 shows a plot of the 
dust temperature as a function of angular distance and 
in units of the inner radius of the dust shell, respectively. 
The inner CDS radius rjn can be derived directly from the 
radial profile, as it is the point at which the dust reaches 
Tin = 1000 K. It corresponds to an angular radius of i^in = 
10.5 mas. 
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Fig. 4. SED (left) and visibility (right) of OH 104.9+2.4 for our best-fitting model. As described in the text, the bolometric 
flux used for the SED is a factor of 3.3 lower than the one used for the SED, in order to take into account the different epochs 
of the ISO and SAO observations. 

Table 2. The derived physical parameters of OH 104.9+2.4 as provided by the best-fitting model. 



Parameter 



Value 



Varied range Impact on SED/visibility 



Effective temperature (black-body) 
Temperature at inner CDS boundary 
Density profile within the CDS 
Relative CDS thickness 
Dust grain distribution function 
Minimum grain size 
Maximum grain size 
Dust types 

Optical depth 



Radius 

Radius of inner CDS boundary 
Radius of / = 10" ''^ /max 
Outer radius of the model 
Mass-loss rate 



Toff = 2500 + 500 K 

rin = 1000 + 200 K 

p{r) (X r-" with n = 2.0 + 0.1 

^ = 10" with p = 5+^ 

MRN, n{a) ex a"' with q = 3.5 

amin = 0.005 + 0.003 ^m 

flmax = 0.2 ± 0.02 fim 

95 % OHM warm sihcates 

5 % DL astronomical silicates 

To.SSfim = 158 ± 7 
T2.2Mm = 6.5 ±0.3 
T9.7;.m = 14.0 + 0.6 

R* ~6OO7?0 =2.79AU 

i?in = 9.1i?* =25.4AU 

Rout° = 5 ■ 10^ i?, = 1.4 ■ 10'' AU 

i?out = 9.1 • 10^ R^ = 2.5 • lO'' AU 

M = 2.18- 10~^ Mq/jt 



1500 . . . 5000 K 

500... 3000 K 

-1...4* 

1...7 

2 ... 5** 

0.0005 ... 0.5 fim 

0.05... 500 ^m 



0.5... 100 



very weak/weak 

moderate/moderate 

moderate/weak 

very Aveak/very wreak (p^3), 
moderate/weak (p^3) 

moderate/strong 
very weak/moderate 
weak/very strong 
very strong/weak 



moderate/strong 



parameter variation includes broken power laws, superwind models and models with radiatively driven winds 
: parameter variation includes test of alternative grain size distributions (see text) 



The bottom right panel of Fig. 5 shows the total optical 
depth as a function of wavelength. As expected, local max- 
ima of T\ are found for the silicate absorption features at 
9.7 ^m (SiO stretching vibrations) and 18 //m (Si02 bend- 
ing vibrations). 

The stellar radius obtained for the best-fitting model 
is "(?* — 1.16mas. This corresponds to i?* ci 600 i?© or 
2.79 AU for the adopted distance of Z) = 2.38 ± 0.24 kpc. 
Therefore, the inner radius of the CDS is i?in = 9.1 i?* 
= 25.4 AU. This value is rather high and results from the 
very high optical depth leading to high dust temperatures, 
so that closer to the central star the dust grains evap- 
orate (see Men'shchikov et al. 2002b). The outer radius 



of the CDS as defined above is at R, 



-10 
out 



1.4 • 10-* AU. From the definition of R~^^ 



i?„„\°/i?in = 550, and since /„ 



''out 



= 5 • 10''' i?* = 

it follows that 

) = 10-10 /„,ax, only 



a fraction of the order of 10^^ of the total flux comes 
from a region with r > i?(7ut • This explains why mod- 
els with Tout /''in > 10'^ do not show much of a change 
with rising Tout /'"in- The value obtained under this as- 
sumption is of the same order as that found by Herman 
& Habing (1985) and therefore seems to give a reliable 
lower limit for the outer diameter of the dust shell. They 
derived Rg = 4.37 ± 0.42 • lO^^ cm = 2.92 ± 0.28 • 10^ AU 
= 6.3 • lO^i?0 ~ 10'^ i?*. However, the modeling as well 
as the observations are both fully allegeable with a dust 
shell that extends from i?in to infinity. 

From the observations, the total luminosity can be de- 
rived according to 



bol 



[10-10 VF/m2] 



3.21 



L 



f D 



[104^0] V[kpc] 



(1) 
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contribution from 
thermal dust emission 
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Fig. 5. Modeling results of OH 104.9+2.4 for our best-fitting model (see Tab. 2). Top left: Fractional contributions to the 
total flux. The contribution from direct stellar light is small in comparison to the contribution of the processes in the CDS. 
Top right: The normalized intensity profile at 2.13 ^m. The sharp central peak corresponds to the central source of radiation. 
Bottom left: Dust temperature as a function of angular distance (lower axis) and of the radius in units of the inner dust shell 
radius (upper axis). The point where Tdust ~ 1000 K indicates the inner radius of the dust shell (^- = 1). Bottom right: 
Wavelength depencence of the total optical depth. 



This leads to a luminosity of Lobs = (1-23 ± 0.12) • 
10"* Lq for tlie bolometric flux at <&iso which is in very 
good agreement with the value Lmod = 1.262 • 10^ Lq 
obtained from our best-fitting model using the Stefan- 
Boltzmann law. 

According to the mass-luminosity relation by Ivezic 
et al. (1995), this corresponds to a central star mass of 
M ~ 1 Mq, which is in line with predictions from stellar 
evolution theory for stars in this evolutionary stage. 

Following Ivezic et al. (1995), the mass-loss rate M can 
be derived from 



Mvc 



L 



TF 



1 - 



1 



(2) 



where 



r = 



-^ grav 

is the ratio of radiation pressure i^rad and gravitational 
force i^grav pei unit volume in the envelope. With rp — 3.2, 
r = 1.7 ■ L from our best-fitting model and the adopted 
outflow velocity Ve = 15km/s from te Lintel Hekkert et 
al. (1991), M = 2.18 • IQ-^ Mg/yr is obtained. 



4.4. Parameter variations 

In this section, we discuss the impact of several input pa- 
rameters on the process of radiative transfer modeling. 
The best-fltting model (see Sect. 4.3) is used as a basis for 
the discussion, and the effect of the variations of the pa- 
rameters in question is outlined. The discussion is divided 
into subsections, in which the change of each parameter 
is discussed separately. To maintain the presentation as 
concise as possible, the variations for minor parameters 
like aniin and the discussion of more complex models (e.g. 
broken power law density distributions, multiple discrete 
dust shells), that did not deliver better fits than the mod- 
els with less assumptions are omitted. 



4.4.1. Dust types 

One of the most important quantities that have to be 
fixed in the modeling process is the chemical composition 
of the dust itself. Here, four different dust types are un- 
der discussion, namely (i) that of the best-fitting model, 
which is mainly the warm, oxygen-deflcient silicates from 
Ossenkopf, Henning & Mathis (1992), (ii) the type given 
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by the same authors as cold, oxygen-rich silicates, (Hi) 
the astronomical silicates by Draine & Lee (1984), which 
are also included as a small fraction in the dust type of 
our best-fitting model (5%), and finally (iv) the a-SiC 
grains tabulated by Pegourie (1988). The results of radia- 
tive transfer modeling obtained for these different grain 
types, with all other parameters fixed to the value of the 
best-fitting model, are presented in Fig. 6a. The models 
with OHM cold silicates give a slightly worse fit to the 
data than the models with the warm silicates. In general, 
the DL silicates provide good model fits, though these are 
not as good as those provided by the OHM warm sih- 
cates. The small fraction of DL silicates was added to the 
dust composition of the best-fitting model to correctly re- 
produce the depth of the 9.7 /xm silicate feature. With all 
other model parameters fixed to the values of the best- 



fitting model (see Sect. 4.3), we need T2.2/ 



2.0 to fit 



the SED for DL silicates, implying excessive flux at the 
shorter wavelengths, but T2.2/^m = 5.0 for a good visibility 
fit. The Pegourie SiC composition seems to be unsuitable, 
as no good overall fit to SED and visibility can be found. 
It is possible to reproduce the depth of the silicate feature 
with an optical depth of T2.2;im = 18.0, but with a great 
vertical offset (which clearly is due to the fact that SiC has 
vibrational modes different from SiO), and the visibility 
can be fitted with T2. 2/^111 = 6.0. 

4.4.2. Effective temperatures 

Figure 6b shows the variations of the model SED and vis- 
ibility with the effective temperature of the central star. 
For a fixed optical depth, a higher value for T^e leads to 
lower flux at both shorter wavelengths and at the 9.7 /im 
silicate feature, and higher flux at the long-wavelength 
tail. The changes are qualitatively similar to those intro- 
duced by a variation of Tjn (see next subsection) but much 
weaker due to the rather large T2.2/jm considered here. The 
main reason for the changes described is the shift of the 
peak-intensity wavelength of the stellar black-body spec- 
trum which induces changes in the radial temperature pro- 
file of the dust: the dust temperature falls of more steeply 
for higher effective temperatures. The effect on the visi- 
bility is rather small in slope and curvature and can be 
compensated by a small change in T2.2^im, as shown for 
Toff = 2000 K in Fig. 6b. The visibilities for the model 
with T2.2 fim = 6.5 and Tcs — 2500 K and the one with 
T2. 2^*111 = 7.0 and Teff = 2000 K are almost identical. 



4.4.3. Temperatures at the inner boundary 

Figure 6c shows the model variations with changing tem- 
perature Tin at the inner dust shell boundary. For r2.2/^m 
as high as in our favored model, decreasing Tin lowers the 
flux at shorter wavelengths, increases the depth of the 
9.7 /im silicate feature, and raises the flux at the long- 
wavelength tail. This behavior is similar to the changes 
produced by increasing r2.2^m- Therefore, changes in the 



SED due to a variation of Tin can, in principle, be compen- 
sated by a corresponding variation of T2.2 /^m, as illustrated 
in Fig. 6c for Tin = 850 K and Tin = 1200 K. On the other 
hand, the depth of the silicate feature constrains the choice 
of the optical depth. For the visibility function, lowering 
Tin leads to a steeper decline at low spatial frequencies 
and a stronger curvature, and therefore a shallower fall- 
off at high spatial frequencies. With an adjusted optical 
depth (as shown for Tin = 850 K and Tin = 1200 K), the 
steepness rises, and the curvature decreases with increas- 
ing Tin. Finally, it can be concluded that even if models 
of similar quality to the best-fitting model can be found 
with different Tin, the value of Tin will be close to 1000 K. 



4.4.4. Reference optical depths 

The optical depth t is the main parameter in this study. 
The variation of this quantity has a strong impact on both 
the SED and the visibility (see Fig. 6d). For the SED, the 
flux at short wavelengths for low optical depths is rela- 
tively high, whereas at long wavelengths it is relatively 
low. In addition, the SiO feature at 9.7 /im is found in 
emission at the lowest values for T2.2/im- For high optical 
depths, more flux is reprocessed to long wavelengths as 
a result of the stronger reddening effects. The most im- 
portant constraint for this parameter is the depth of the 
9.7 /im absorption feature, which is too deep for very high 
optical depths. The visibility mostly keeps its curvature, 
but steepens for higher optical depths due to the stronger 
obscuration of the central source. The steepening shows a 
roughly linear dependence on the optical depth in the do- 
main under consideration. The best value for a simultane- 
ous SED and visibility model is found to be T2.2/im = 6.5. 



4.4.5. Grain sizes 

In Fig. 7a, models assuming an MRN grain-size distribu- 
tion with different maximum grain size amax are shown. 
All other parameters are fixed to the values of the best- 
fitting model presented in Sect. 4.3. Concerning the SED, 
flmax provides a strong constraint for the scattering and 
absorption efficiency. The averaged (over the size distri- 
bution) scattering efficiency per unit volume of the grains 
QsQe-jy oc c? steeply declines with increasing wavelength, 
and can be neglected above a certain wavelength, depend- 
ing on flniax (for an example, we refer to Fig. 5, top left 
panel of the best-fitting model). The size-averaged ab- 
sorption efficiency depends on the grain size only at short 
wavelengths and becomes independent of amax when A 3> 
o (see Fig. 8). Therefore, the fluxes observed at the short- 
est wavelengths place an upper limit on the maximum 
grain size. The depth of the 9.7 /im silicate absorption fea- 
ture itself leads to the prediction that Cmax < 0.50 /im for 
the current configuration. In the long wavelength regime 
(A > 100 /im), it is possible to find a model that fits bet- 
ter than the best-fitting model, but the improvement of 
the fit quality is relatively small. The slightly better fit 
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Fig. 6. SED (left) and A"-band visibility (right) of OH 104.9+2.4 for the best-fitting model (thick solid hne, 95 % OHM warm 
silicates, 5 % DL silicates) and the impact of several parameter variations, a: Variation of the dust composition. Apart from 
the best-fitting model models with other dust types and values of r2.2,jm (see labels) are shown. All other model parameters 
are fixed to the values of the best-fitting model. The models for DL silicates with T2.2fim = 2.0 and r2.2fim ~ 5.0 are given as 
examples for models that fit the SED and visibility with this specific dust composition. The same comparison is shown for the 
Pegourie (1988) SiC grains with r2.2^im = 6.0 and T2.2 in-n = 18.0, respectively, b: Variation of Teff. The best-fitting model is the 
one with T^s = 2500 K. For comparison, a model with Tcff = 2000 K and T2.2,jm = 7.0 is shown, to illustrate that changing the 
optical depth can make up for the variations caused by Toff (see text for further details), c: Variation of Tin. The best-fitting 
model is labeled with Tin = 1000 K. For comparison, models with Toff = 850 K and r2.2fim ~ 4.5 as well as Toff — 1200 K and 
T2.2/jm = 8.0 are shown, to illustrate that a change of the optical depth can also make up for the variations caused by Tin. d: 
Variation of r2.2nm. The best-fitting model is labeled with r2.2;jm = 6.5. 
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Fig. 7. SED (left) and A'-band visibility (right) of OH 104.9+2.4 for the best-fitting model (thick solid line, 95 % OHM warm 
silicates, 5% DL silicates) and the impact of several parameter variations, a: Variation of the maximum grain size flmax- The 
best-fitting model is the one with amax = 0.20 /.im. b: Variation of the grain size distribution. The best-fitting model is the 
one with the MRN distribution. For comparison, models with the KMH distribution and r2.2fim = 2.0 and r2.2fim ~ 9.0 are 
plotted as examples of good visibility and SED fits, respectively, c: Variation of the dust shell thickness. The best-fitting model 
is the one with ^-^ — 10^. d: Variation of the slope of the radial density distribution. The best-fitting model is the one with 
p{r) oc r~^°. In addition to the p oc r~" distributions a distribution based on a radiatively driven wind model is shown (see 
text for details). 
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der examination. In this study, variations of the standard 
MRN function have been tested. 

Apart from the MRN distribution, the distribution 
from Kim, Martin, & Hendry (1994, hereafter KMH) is fre- 



quently used. It assumes n{a) oc a 



-a/ciQ 



X[\S!m] 

Fig. 8. Total optical depth spectrum for different maximum 
grain sizes (see labels). 

might be due to a fixed _Fboi in combination with the fact 
that the models in question become optically thin at short 
wavelengths (Fig. 8), since the large grains engross large 
fractions of the dust mass. In this case, a large fraction of 
i^boi is emitted at short wavelengths, resulting in lower flux 
at long wavelengths. In addition, the problems mentioned 
in Sect. 4.3.1 may play a role. 

The K'-hand visibility is very sensitive against scat- 
tering and thus depends strongly on the assumed grain 
sizes. Increasing flmax has a similar effect to increasing 
the optical depth and results in a stronger decline of the 
visibility function. As an MRN-like distribution declines 
very steeply towards larger grain sizes, most of the par- 
ticles in the model remain small when Omax is increased. 
But, due to the strong dependence of the scattering and 
absorption properties on the grain sizes at 2.12 //m, the 
impact of these few extra large grains on the visibility 
is comparable to that of the large bulk of small grains 
(Kriiger & Sedlmayr 1997). The curvature of the visibility 
changes its sign somewhere around a^^x < 0.20 /im. This 
behavior reflects the change of the spatial intensity dis- 
tribution. At large radial offsets 6 from the central star, 
7(6) oc h^^ holds, since the optical depth along the line of 
sight becomes small (Jura & Jacoby 1986). 

Regarding the wavelength dependence of the total op- 
tical depth, the choice of Omax is only of minor impor- 
tance beyond a certain range of values (see Fig. 8). For 
Omax < 1 M^i and wavelengths larger than 5 /^m, the im- 
pact of the grain size is limited. The shape of the function 
is almost the same, and the offset can be fixed by adjusting 
the optical depth of the models. If the grains become too 
large and too few, the function approximately converges 
towards a power law, since the characteristics of the grains 
then play a minor role. 

4.4.6. Grain-size distributions 

Instead of the standard MRN grain-size distribution, sev- 
eral other GSD functions can be adopted for the dust un- 



with q = 3.5, 



for a > flfninj where Oq has to be fixed for the exponential 
falloff, which replaces the sharp upper cutoff of the MRN 
distribution. The value chosen here is oq = 0.17 /im, which 
should lead to a distribution comparable to the standard 
MRN. The results are shown in Fig. 7b. For T2.2tim = 6.5, 
the fit is not good. For a good SED fit, r2.2^m = 9.0 is 
needed. However, we need r2.2/im — 2.0 for a good visibil- 
ity fit. Therefore, models with a KMH grain-size distribu- 
tion provide a good fit of either the SED or the visibility, 
but not both simultaneously. 

In addition to KMH we implemented in the DUSTY 
code the GSD from Witt et al. (2001). Here, the MRN 
function remains unchanged, but a certain break point 
^brcak is introduced, where the slope of the exponen- 
tial decline q changes. We first used the original values 
qi = 3.5, 92 = 4.0, flniin = 0.005 ^m, abrcak = 0.5 ^m and 
a-max = 2.0 /xm given by Witt et al. (2001). Obviously, re- 
sults obtained with this parameter set are not directly 
comparable to the best-fitting model, due to different 
grain sizes. With abrcak = 0.15 ^m and amax = 0.22 /im, 
an SED fit of comparable quality to our best-fitting model 
is obtained, while the visibility fit is worse at intermedi- 
ate spatial frequencies. Since reasonable improvements all 
converge towards the parameters of our best-fitting model, 
our assumption about the MRN GSD is justified. 

4.4.7. Dust-shell thickness 

In Fig. 7c different relative dust-shell thicknesses are ex- 
amined. Except for the model with ^^'^^ = 10, the SED 
fits are almost the same up to A ~ 30 /xm, while the model 
visibility remains essentially unchanged. Taking into ac- 
count the poor visibility fit, dust shells with ^-^^^ <C 10^ 
can be ruled out. Above 30 /im, the SED falls off more shal- 
lowly for rising geometrical thickness, providing a better 
fit to the observational data. The additional cold dust in- 
troduced by the larger outer boundary increases the far 
infrared flux but does not affect the SED at shorter wave- 
lengths. The only difference between the models ^^^^^ = 
lO'' and ^^^^^ — 10^ is a small deviation for A > 120 um. 
Thus, from our analysis the model with ^^'^^ = 10^ is 
slightly favored, whereas the relative dust-shell thickness 
turned out to be a non-crucial parameter in the modeling 
process (see also discussion above about R^ul^)- 



4.4.8. Density distributions 

So far, it was assumed that p oc r~^ holds within the 
dust shell. Assuming the outflow velocity Vc is the same 
as for the standard model and constant within the CDS 
and a constant dust-to-gas ratio, any density slope shal- 
lower/steeper than p oc r~^ means that the stellar mass- 
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loss rate M is non-stationary and decreases/increases with 
time. Figure 7d shows the result of variations in the den- 
sity distribution. Since the overall agreement between ob- 
servations and models worsens for p ex r~^ with a; 7^ 2, 
our standard density profile is justified. 

For AGB stars, the envelope expansion is driven by 
radiation pressure on the dust grains. Therefore, a sim- 
ple dust density profile p{r) oc r^^ may not be the most 
realistic one. If the variation of the flux-averaged opac- 
ity with radial distance as well as the grain drift is neg- 
ligible, an analytic solution for the problem of a radia- 
tively driven wind can be found (Ivczic et al. 1999) which 
might be more appropiate for modeling cool evolved stars 
such as OH 104.9-1-2.4. In Fig. 7d such a model is shown, 
whose other parameters are identical to those of our best- 
fitting model. The SED fit agrees with the observations at 
most of the wavelengths discussed. Only for short wave- 
lengths (A < 3 /im) and long wavelengths (A > 80 pxn) is 
there a noticeable deviation from the observational data. 
Unfortunately, no better model fit of the long wavelength 
regime could be obtained compared to our best-fitting 
model without a radiatively driven wind. In terms of the 
visibility, some models provide fits of comparable quality 
to our best-fitting model shown in Fig. 7d or even slightly 
better ones (with a worse SED fit). After all, both model- 
ing approaches provide models of similar quality with oth- 
erwise identical parameters, while our best-fitting model 
is still slightly favoured. 

5. Summary and conclusions 

We present diffraction- limited 2.12 /im speckle interfcro- 
metric observations of the circumstcllar dust shell around 
the highly obscured type II-A OH/IR star OH 104.9+2.4. 
The resolution achieved with the SAO 6 m telescope is 
74 mas, which is sufficient to fully resolve the CDS at this 
wavelength. From an azimuthally averaged 1-dimcnsional 
Gaussian fit of the visibility function, the diameter was de- 
termined to be 47 ± 3 mas (FWHM), which corresponds 
to a diameter of 112 ± 13 AU for an adopted distance of 
D — 2.38±0.24kpc. The reconstructed 2-dimcnsional vis- 
ibility shows no major deviation from spherical symmetry. 
In order to investigate the structure and the proper- 
ties of the CDS, we performed 1-dimensional radiative 
transfer calculations to model simultaneously the spectral 
energy distribution (ISO SWS/LWS spectrum), as well 
as the 2.12 /im visibility function from our observations. 
The obtained dust shell properties correspond to a stan- 
dard single-shell model (uniform outflow). Multiple com- 
ponent models do not match the observations better as the 
simple uniform outflow models, while they are certainly 
not excluded with the remaining uncertainities. Among 
the ^ 10^ models calculated, the ISO SED, which ap- 
pears to correspond to the minimum pulsation phase of 
OH 104.9+2.4, was best reproduced by a model where the 
inner rim of the dust shell (i.e. the dust condensation zone) 
has a dust temperature of Tjn = 1000 K and an angular di- 
ameter of 10.5 mas corresponding to 9.1 i?* (Fig. 5, bottom 



left panel). At 2.12/im, this zone is not limb-brightened, 
since the model is optically thick (r2.2^m = 6.5). 

During the entire model grid calculations, the validity 
of the parameters given above as well as of other param- 
eters such as dust temperatures, grain sizes, and effec- 
tive temperatures was comprehensively investigated. The 
grain-size distribution was found to agree with a standard 
MRN distribution function, n{a) ex a^^'^, with a ranging 
between Omin = 0.005 /im and Cmax = 0.2 /im, which shows 
that the influence of more complex processes like the dust 
formation process itself on the model seem to be rather 
little. A two-component dust composition consisting of 
95 % OHM warm silicates and 5 % DL silicates was found 
to match the observations best. Assuming a single black- 
body spectrum for the spectral shape of the central source, 
it is best represented by a model with Tgff = 2500 K. 
However, moderate variations of Toff only slightly affect 
the model results. From the best fitting model, we derived 
a bolometric flux at Earth of Fboi — 6.99 • 10~^^ W/m^, in 
accordance with the ISO observational data for the SED 
at ^iso = 0.5 = $min. This bolometric flux corresponds 
to a luminosity of L = (1.23 ± 0.12) • 10^ Lq at a distance 
of I? = 2.38 ± 0.24 kpc. The radius of the central star is 
found to hedi, — 1.16 mas or i? ~ 600 Rq, and the central- 
star mass was derived to be M 2± 1 Mq using the method 
of Ivezic et al. (1995). Though the density distribution 
was varied greatly, the general p{r) oc 1/r^ dependence 
was proven to best flt the observations. The present study 
matches the observations for a broad wavelength range 
(from 1 /im to 220 /im). 

As regards our visibility at 2.12/im, the difference 
in pulsation phase between the SAO observations and 
the ISO measurements had to be taken into account. 
According to the observations, the bolometric flux be- 
tween the 1996 ISO SED phase and the 2002 SAO visibil- 
ity phase differs by a factor of 3.3. Therefore, this factor 
had to be applied to the bolometric flux of the visibil- 
ity model in order to simultaneously fit the SED and the 
visibility. 

As concerns the study of the parameter space, the 
grain size turned out to be the key parameter in adjusting 
the visibility fit while preserving the quality of the SED fit. 
Both the slope and the curvature of the visibility depend 
sensitively on the assumed grain radii. The dust mass-loss 
rate M is well constrained by the shape of the SED at 
longer wavelengths and especially by the shape of the sil- 
icate absorption feature. For given optical constants, the 
value of M, as derived from the match of the feature, is 
not very sensitive to changes in the input parameters. The 
derived value of M = 2.18 • 10^^ Mq/jx is of the same or- 
der of magnitude as the value of M = 5.58 • lO~^M0/yr 
given by Heske et al. (1990). The shape of the observed 
visibility and the strength of the silicate feature constrain 
the possible grain radii and optical depths of the model, 
which are found to be the same values that provide the 
best fit to the SED. 

While our model of OH 104.9+2.4 simultaneously fits 
the observed SED and visibility very well, the modeling 
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included some simplifications. The central source is not 
a perfect black body, as the dynamical stellar atmosphere 
exhibits a complex structure and heterogeneous chemistry. 
In addition, the feature at 18 /im could not be reproduced 
very well in the SED, due to the simplified dust compo- 
sition, and the wavelengths above 100 /im could only be 
fitted to a certain extent. This is probably due to the 
fact that DUSTY allows for the entry of a grain-size dis- 
tribution, but the code averages dust properties over the 
entire size distribution before performing any further cal- 
culations (see e.g. Wolf 2003). In our case, this leads in 
particular to an absence of the very large grains found in 
the tail of the grain-size distribution, which have very high 
and almost constant absorption and scattering efficiencies. 
By choosing larger grain sizes, we obtain a much better 
SED fit for A > 50 /im, as the emission from the large 
grains, which have a flatter absorption efficiency and a 
higher scattering efficiency as smaller grains, then traces 
the shape of the SED at the long wavelengths, where the 
dust shell is otherwise optically thin. Unfortunately, this 
also leads to an optically thin solution for shorter wave- 
lengths due to a lack of small grains, as most of the dust 
mass is in this case utilized by large grains. 

For future considerations, it would be very useful to 
have additional K-hand photometry to secure the derived 
value for the pulsation phase, as well as photometry in 
other bands to get better insight into the wavelength de- 
pendence of the pulsation amplitude. More visibility mea- 
surements, e.g. in the J, H, N, and V bands, would pro- 
vide additional information on the structure of the dust 
shell, while additional visibilities in the K' band at an- 
other pulsation phase would allow for more detailed con- 
clusions about the response of the dust shell to the vari- 
ation in the central object, as there may be light-travel 
time effects involved due to the large scale of the CDS 
itself. 
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